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ABSTRACT 



Context. We report the detection of five high-velocity stars in the core of the globular cluster NGC 2808. The stars lie on the red giant 
' branch and show total velocities between 40 and 45 km s~' . For a core velocity dispersion <r c = 13.4 km s~' , this corresponds to up to 

3.4o" c . These velocities are close to the estimated escape velocity (~ 50 kms~' ) and suggest an ejection from the core. Two of these 
stars have been confirmed in our recent integral field spectroscopy data and we will discuss them in more detail here. These two red 
(-H . giants are located at a projected distance of ~ 0.3 pc from the center. According to their positions on the color magnitude diagram, 

both stars are cluster members. 

Aims. We investigate several possible origins for the high velocities of the stars and conceivable ejection mechanisms. Since the 
velocities are close to the escape velocity, it is not obvious whether the stars are bound or unbound to the cluster. We therefore 
consider both cases in our analysis. 
Methods. We perform numerical simulations of three-body dynamical encounters between binaries and single stars and compare the 
resulting velocity distributions of escapers with the velocities of our stars. If the stars are bound, the encounters must have taken place 
when the stars were still on the main sequence. We compare the predictions for a single dynamical encounter with a compact object 
with those of a sequence of two-body encounters due to relaxation. If the stars are unbound, the encounter must have taken place 
recently, when the stars were already in the giant phase. 

Results. After including binary fractions and black-hole retention fractions, projection effects, and detection probabilities from Monte- 
Carlo simulations, we estimate the expected numbers of detections for all the different scenarios. Based on these numbers, we conclude 
that the most likely scenario is that the stars are bound and were accelerated by a single encounter between a binary of main-sequence 
stars and a ~ 10 M Q black hole. Finally, we discuss the origin of previously discovered fast stars in globular clusters, and conclude 
that the case of NGC 2808 is most likely a representative case for most other detections of fast stars in globular clusters. We show 
that with the present analysis we are able to explain high-velocity stars in the clusters M3 and 47 Tucanae with simple dynamical 
f^*) . encounters. 

Key words, globular cluster: individual (NGC 2808) - stars: kinematics and dynamics - stars: high-velocity 
1 . Introduction being field stars is rather low. iMevlan et al.l(ll99ll) found a simi- 
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!_| , lar case for the globular cluster 47 Tucanae. The stars which they 

03 . High-velocity stars in globular clusters have been a puzzle found haye radial velocities of _ 36 7 ^ s -i and 32A ^ s -i in 

to astronomers since the beginning of velocity measurements. ±e reference frame of the cluster which correS p 0n ds to 4.0 and 

Especially finding fast giant .stars in the^ore of globular clusters 3 6 timefj the CQre vdoci di rsion = 9 .1 km s 1 . Despite 

has drawn attention. |Gunn & Gnffin| dl979j) found two fast stars ±e fact ^ ^ 1qw emic yeloci of ±e dugter doefj nQt ^ 

which they called 'interlopers in the globular cluster M3. The 1qw for an accurate kinematic statemen t on the membership of 

stars are located in the core of the cluster, both at a projected these tWQ starSi their Won fa ±e co i or . magnitude diagram 

distance of about 20 from the center. With radial velocities of and ±e M h Galactic of 4? Tucanae (b = _^ 0) M 

17.0 km s- 1 and -22.9 km s" 1 relative to the cluster, they move for membership 

with 3.5 and 4.5 times the central velocity dispersion (cr c = 4.9 

kms 1 ), respectively. Due to the high systemic velocity of the Plausible mechanisms to explain these high-velocity stars 

cluster itself (V r ~ -147 kms" 1 ) the possibility of these stars are ejection from the core by recoil from an encounter be- 

tween a single star and a binary, between two binary stars, be- 

* Based on observations collected at the European Organization for tween a binary and an intermediate-mass black hole (IMBH) at 

Astronomical Research in the Southern Hemisphere, Chile (083.D- the center or even from encounters of stars with an IMBH bi- 

0444). nary (e.g. iGualandris eTaIll2004l l2005t iBaumgardt et al.ll2006t 
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Fig. 1. ARGUS field of view (left) and velocity map of NGC 2808 (right). Visible as bright blue spots are the two high-velocity stars 
FS1 and FS2. Marked with the dashed blue line is FS4, the star which does not dominate its spaxel and therefore is not detectable. 
The blue cross marks the center of the cluster. 
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Fig. 2. Kinematic properties of the five high-velocity stars. Left: The histogram of the Fabry-Perot data in the core of NGC 2808. 
Marked with colors are five high velocity stars FS1 to FS5. Right: The escape-velocity profile (solid line) and its uncertainties 
(dotted line) obtained from the density profile overplotted by the total velocities and projected positions of the five fast stars. 



iGvaramadze et al.1 12009: Mapelli et al.1 120051) . It has been re- 
cently suggested that globular clusters may c ontain central black 
holes in the mass range of 10 3 - 10 4 M G dMiller & Ham ilton 
l2002t iBaumgardt et al.ll2005l) which fall on top of the relation 
between the velocity dispersion and black hole mass seen fo r 
galaxies (e.g. iFerrarese & Merrittl kOOOl iGebhardt etaT1l2000h . 
Kinematic signatures for intermediate-mass blac k holes have 
been found in the globu lar clusters cj Centauri (jNovola et al.l 
2008, [20loh. Gl in M31 (IGebhardt et al.l l2002l |2005|) and NGC 
6388 dLiitzgendorf et al.ll2.Ql ll) . High velocity stars are predicted 
in globular clusters which host intermediate-mass black holes. 
Therefore it is important to study such stars in more detail when 
observed in a globular cluster. 

The globular cluster NGC 2808 is of high interest in many re- 
gards. From the photometric s ide, it has a complex extended hor- 
izontal branch (iHarrisI ll 974l iFerraro et al.l 119901) which shows 
puzzling discontinu ities in the s tellar distribution along its blue 
tail. Furthermore, iPiotto et aD (120071) found a triple main se- 



quence after accurate photometric and proper motion analy- 
sis with deep Hubble Space Telescope (HST) data. This in- 
dicates the existence of three sub-populations with an age of 
~ 12.5 Gyr and different metallicities. The kinematic properties 
of NGC 2808 are also unusual. By analyzing Fabry-Perot data, 
Geb hardt et al.l d2012l) found rotation in the outer parts of the 
cluster with a maximum rotation velocity of V lot ~ 5 kms -1 . 
The latest radial velocity measurements of the cluster core by 
iLiitzgendorf et al.l d2012l) have shown that NGC 2808 most prob- 
ably does not host an intermediate-mass black hole: the best fit 
model is consistent with no IMBH and gives an upper limit of 
M fiH <lxlO 4 M . 

In this work we report on the discovery of five high-velocity 
stars in NGC 2808 and discuss several explanations and ejec- 
tion scenarios for the two stars found in our integral field spec- 
troscopy observations. In Section [2] we explain the observations 
and the method of obtaining single star velocities with integral 
field spectroscopy. Further, we introduce a second data set used 
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for velocity measurements taken from a Fabry-Perot instrument. 
At the end of the first Section we determine the escape velocity 
of NGC 2808. Section[3]is dedicated to the discussion of alterna- 
tive explanations for the observed high velocities including fore- 
ground stars, binaries and atmospheric active stars. Section|4]de- 
scribes the analysis of the Maxwellian velocity distribution and 
the probabilities of these stars resulting from relaxation while 
Section[5]describes the dynamical three-body simulations in or- 
der to find the most likely acceleration scenario for the two fast 
stars. In Section[6]we briefly discuss other cases of fast star de- 
tections in globular clusters. Finally, we summarize our results 
and list our conclusions in Section|7] 



2. Radial velocities 

The data we use were obtained with the GIRAFFE spec- 
trograph of the FLAMES (Fiber Large Array Multi Element 
Spectrograph) instrument at the Very Large Telescope (VLT) in 
ARGUS (Large Integral Field Unit). F or more details on the ob - 
servations and reductions we refer to ILiitzgendorf et aD (t2012h . 
From these data we construct a velocity map which is shown 
in Figure Q] together with the corresponding field of view from 
the HST image. The two fast stars are visible as bright "blue" 
spots and are labeled with FS1 and FS2 in the maps. To estimate 
their velo cities, we use the output of a shot noise routine de- 
scribed in ILiitzgendorf et aTl (120121) . This routine computes for 
every spaxel the number of stars that contribute to the light of 
the spaxel and indicates which spaxels are dominated by a single 
star. With this information we find the spaxel to which our two 
(bright) fast stars contribute with more than 80% of the light and 
use the radial velocities derived for these spaxels as the ones of 
the stars. These velocities are co rrected for the systemic v elocity 
of the cluster (V r = 93.6 km s' 1 ILiitzgendorf et al.ll2012l) . 

A sec ond measurement is o btained from the dataset ob- 
served by Gebhar dt et al.l (1201 2|) with the Rutgers Fabry Perot 
on the Blanco 4-m telescope at Cerro Tololo Inter-American 
Observatory (CTIO). This dataset contains 3634 stars and covers 
the center of the cluster as well as regions up to 4 arcmin radius. 
We match this dataset with our photometric catalog from the 
HST image and identify the two fast stars. Despite crowding in 
the inner regions, the two stars are resolved and can be identified 
in the Fabry-Perot dataset. In addition to FS 1 and FS2, we find 
three more stars with velocities values between 36.0 kms -1 and 
42.0 km s _1 inside the core radius labeled with FS3, FS4 and 
FS5. Table Q] lists the velocity measurements of all the stars as 
well as other properties derived from the HST data. Besides FS 1 
and FS2 only FS4 is also located inside the field of view of the 
combined ARGUS pointing (see Figure [TJ. However, due to its 
low contribution of ~ 55 % to the single spaxel we measure a 
much lower velocity in the IFU data than with the Fabry-Perot. 
The remaining stars lie outside the IFU pointing. Figure |2]shows 
the histogram of the radial velocities inside the core from the 
Fabry-Perot data set. The five fast stars are color-coded and vis- 
ible as outliers of the distribution. 

In addition to the spectroscopic data sets we also received 
proper motions for four out of the five stars (private communi- 
cation, Jay Anderson). As listed in Table Q] the proper motions 
indicate that all stars are cluster members and increase the total 
velocities up to 40 - 45 km s _I . FS3, unfortunately, lies outside 
the proper motion dataset and therefore can only be assigned a 
lower limit of velocity. Since our measurements are only com- 
plete for FS 1 and FS2 we decided to limit the following analysis 
to these two stars. 



Figure [3] shows the photometric data points of the Fabry- 
Perot for FS 1 and the best fit to its spectrum. The data of FS 1 is 
not well reproduced by the model (solid line). In order to com- 
pare velocities from both datasets the velocities are corrected 
for the systemic motion of the cluster by computing the aver- 
age velocity of the central region. The results from both data 
sets show high velocities for these stars. However, for FS1 the 
two measurements differ by 7 kms -1 , i.e. more than lcr, from 
each other. The star is close to the saturation limit and therefore 
causes difficulties when measuring the velocity. Also, measuring 
individual velocities with integrated light can lead to larger er- 
rors since other stars always contribute some amount of light to 
the analyzed spectrum. For FS2 the measurements agree within 
the error bars. Thus, the difference for FS1 is probably not due 
to a systematic shift. The discrepancy for FS 1 can have several 
reasons: a) the velocity of this star is variable. This suggests a 
binary or long period variable and is discussed in more detail in 
Section [3] b) At least one of the measurements is affected by a 
large error. Considering the quality of the fit of the Fabry-Perot 
data we consider the spectroscopic measurement more accurate. 

As a next step, we calculate the escape velocity of the cluster 
as a function of radius. This value gives the information needed 
to decide whether the stars are bound or unbound. We obtain an 
escape velocity profile by parametrizing our surface-brig htness 
profil e with a multi-gaussian expansion (MGE, Emsellem et al.l 
1994) as implement ed in the Jeans Anisotropic Models (JAM) 
bv lCappellaril (12008 ) . The parametrization allows an easy way of 
deprojecting the profile after multi plying it with the derived M/L 
profile from our kinematic fits (see Liitzgen dorf et al.l2012l) . The 
gravitational potential generated by the deprojected de nsity of 
stars is given by Equation (39) of lEmsellem et al.l d!994l) : 

, — r A MfHj{u) 
<f>(R,z) = -tJ2/ttG I 2j du > (!) 

o 7=1 1 

where G is the gravitational constant, Mj and crj the enclosed 
mass and the dispersion of each of the N Gaussians, respectively 
and with 

<Hj(u) = 1 1 V A (2) 

^1~(1~^) M 2 

a function of the integration variable u, the deprojected axial ra- 
tio < qj < 1 ., and the cylindrical coordinates R and z. 

The escape-velocity profile is then evaluated by calculating 
the difference in the potential at a radius o f r and the potential at 
the tidal radius (r, = 35 pc, lHarri si [19961) as given by the equa- 
tion: 

V esc (r) = V2 (%) - <D(r)) (3) 

The resulting profile is shown in Figure Overplotted are 
the velocities of the five stars at their projected radii. Since the 
radii are projected, they only give a lower limit on the actual dis- 
tance from the center. Considering the uncertainties of the stellar 
velocities and of the escape velocity from the cluster, we find that 
star FS4 is consistent with being both bound and unbound. The 
remaining stars appear to be bound,but given the large velocity 
uncertainties we will also investigate ejection mechanisms in the 
unbound regime. 
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Table 1. Radial velocities and photometric properties of the two fast stars FS1 and FS2. The table lists the radial velocity mea- 
surements from the IFU dataset (Vifu) an d the velocities measured from the Fabry-Perot data (Vfp)- niv and nt] are the V and I 
magnitudes of the stars and r cen is the projected distance from the cluster center. 
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Fig. 3. Fabry-Perot photometric points for FS 1 . Overplotted is 
the best fit of the model for the Ha line. 



3. Possible explanations 

Despite the uncertainties in the exact velocity, both stars lie 
more than 3<x e above the sy stemic velocity of the clu ster, where 
cr c = (13.4 + 0.2) kms- 1 (Liitz gendorf et al.l 120121) represents 
the central velocity dispersion. Before investigating dynamical 
acceleration scenarios for the stars, we consider possible alter- 
native explanations for the high velocities. 

3.1. Foreground stars 

The first possibility we consider is that the two stars are field 
stars and therefore moving with a different velocity than the clus- 
ter. We identify the two stars in our HST image (see Figure [1] 
and study their positio n on the color-magnitude diagram (CMD) 
which was obtained in Liitzg endorf et al.l(l2012l) . Figure|4]shows 
the position of all the stars in the CMD. Both stars are located 
on the giant branch, which suggests cluster membership. The 
three other stars are also most likely cluster members. FS1 is 
the brightest star in our pointing and sits at the top of the giant 
branch. The fact that the brightest star in our dataset is also the 
fastest is suspicious. This star is most likely a long period vari- 
able star and will be discussed in more detail in Section l3~3l The 
second star (FS2), however, lies in the central part of the giant 




Fig. 4. Color-magnitude diagram of NGC 2808 overplotted with 
the positions of the high-velocity stars FS 1 to FS4. All of them 
are located on the red giant branch or the horizontal branch and 
therefore most likely cluster members. 



branch about one magnitude fainter than FS 1 and does not seem 
to show any peculiar photometric properties. 

The fact that both stars lie on the giant branch of the clus- 
ter does not strictly exclude them from being foreground stars. 
Stars at a distance of a few hundred pc may also appear super- 
posed on the giant branch of NGC 2808. In the case of FS2, a 
lower main sequence foreground star would be projected up on 
the giant branch if its apparent magnitude was increased by 8 to 
9 magnitudes. This would imply that the star is located 40 - 60 
times closer than a giant of similar apparent magnitude. That is, 
at a distance of 240 - 150 pc from us. FS 1 is even more extreme. 
With a color (B — V) ~ 1 .9, an equivalent dwarf star would have 
an apparent brightness increased by more than 12 magnitudes. 
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This would mean that the star is 240 times closer to the sun than 
its giant equivalent and therefore at 40 pc distance from us. 

In general the probability of contamination through galac- 
tic field stars is rather low. A ccording to a galactic model by 
iRatna tunga & Bahcall (1985), the estimated number of field 
stars per square arcminute projected onto the central regions of 
NGC 2808 is 3.4 x 10~ 2 for 0.8 < (B - V) < 1.3 (FS2) and 
1.8 x 10 3 for (B - V) > 1.3 (FS1) The stellar popula tion syn- 
thesis model of the Milky Way bv lRobin et alj d2004l) can also 
be used to predict the number of foreground stars expected in 
this part of the sky. For a field of view of 0.5 arcmin 2 at the po- 
sition of I = 282.19°, b = -11.25°, the model does not predict 
any stars in the observed magnitude ranges (my = 13.5 - 14.5 
mag, mi = 11.5 - 13.0 mag). Also, the high radial velocity of 
NGC 2808 makes it easier to distinguish between field stars and 
cluster members. 

Finally, neither star shows peculiar motions in the plane of 
the sky and the tangential velocities fall inside the velocity dis- 
tribution of the cluster (see Table [TJ. This strongly argues for 
membership and a dynamical connection between the stars and 
cluster. Considering all these facts, we conclude that both stars 
are giants and associated to the cluster. 

3.2. Binaries 

The binary fraction in typical globula r clusters is about 
10 % (e.g. [ Rubenstei n~& Bailvnl fl997t iDavis et alJ 120081: 
iDalessandro et al.ll201 ll) . Due to their higher dynamical mass, 
binaries tend to sink to the cluster center and become more con- 
centrated in that region. This is why it is likely that binaries are 
present in our field of view. The possibility that these two fast 
stars are binaries, however, is rather low. 

The fact that we measure very similar radial velocities at dif- 
ferent epochs (IFU and Fabry-Perot datasets) with a time dif- 
ference of 15 years already speaks against a binary scenario. 
Within this time period we should have seen large differences 
between the two measurements if the velocities were periodi- 
cal. If we assume that the radial velocities that we measure for 
our two stars are their orbital velocities, a ~ 1 M Q companion 
would imply a separation of ~ 0.5 AU. The two stars are located 
on the giant branch and their radii are of about 7 R (FS 1 ) and 
70 R (FS2) (most likely lower limits), which are comparable to 
the presumed separation. This would indicate that these systems 
are in deep contact and that we should see signatures of accre- 
tion and mass transf er in the form of X-ray and UV excesses. 
IServillat et al.l (2008) obtained X-ray observations with Chandra 
for the core of NGC 2808. Comparing their observations with 
our image, we see no match between our two stars and any of the 
X-ray sources detected. We also investigate a near and a far UV 
image obtained with the Space Telescope Imaging Spectrograph 
(STIS) on board the Hubble Space Telescope (Program ID: 
8511). The stars are visible in the near UV but no source was 
detected in the far UV image. We conclude that these stars show 
no excess in either X-ray and UV which excludes almost all sce- 
nar ios of contact binaries. 

ISommariva et al.l ( l2009t) investigate the binary fraction of the 
globular cluster M4 via spectroscopy. They find 57 binary star 
candidates out of 2469 observed stars with 4 candidates inside 
the core radius. Looking at their table of candidates we find 
7 binaries which show velocities > 20 kms -1 in one of their 
orbital phases. Since the velocity dispersion of NGC 2808 is 
over twice as high as the one of M4 (~ 4 kms -1 ), total ve- 
locities around 40 kms -1 would be possible if the binary itself 
approached with a velocity > 20 kms" 1 . However, all high ve- 
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Fig. 5. Spectra of the two fast stars extracted from the ARGUS 
integral field unit in the Calcium Triplet region. The lines of both 
stars are narrow and confirm the assumption that both stars are 
giants. 



locity binaries in M4 lie in the upper main sequence and the 
turn-off region. There is no high velocity binary observed on the 
giant branch. To make a quantitative statement, we estimated the 
expected fraction of binaries with these properties and velocities 
higher than 45 km s _1 with binary population synthesis. 

We assume a binary fraction of 10 % and use a log-normal 
orbital period distribution with a mean of (log (P/days)) = 4. 8 
and a sigma of 2.3 as obtained by Dua uennov & Mavoi (Il99ll) . 
The binary components are chosen randomly from an evolved 
(over 12 Gyr) Kroupa initial mass function (IMF) to account for 
the age of NGC 2808. For the eccentricities we assume a ther- 
mal distribution of /(e) = 2e. Thus, for every binary we com- 
pute M\,Mi,P and e. From P we derive the semi-major axis 
a and place both components (M\,Mi) at random positions in 
their orbits. We then distribute the binary orientation randomly 
in space and add systemic velocities drawn from a velocity dis- 
tribution with <x = 13.4 kms~' . Finally we observe the binaries 
from a random position and count how often a giant star or main- 
sequence star with M > 0.8 M has velocities |V| > 45 km s _1 . 
As additional boundary conditions we required that the pericen- 
ter distance of the binary is larger than the sum of the stellar 
radii, i.e. a{\e) > (R*i + R*2) 

Among the 1000 bright stars in our pointing we expect 
0.0054 stars like FS2 (R* ~ 7R ) and only 0.00012 stars like 
FS 1 {R* ~ 70 R Q ) if their velocities originate from orbital bi- 
nary velocities. To summarize, it is possible but highly unlikely 
that the high velocities we observe are due to orbital velocities 
in binary systems. 



3.3. Atmospheric active stars 

Another possibility is that these stars have active atmospheres, 
such as strong stellar winds or expanding shells which we 
would see as approaching velocities. According to its loca- 
tion in the CMD, FS1 is most likely a Mira star. Studies have 
shown that the amplitudes of Mira pulsations do not exceed 30 
kms" 1 (e.g.lHinklelll978HHinkle et al Jl 1982b [Hinkle et alJll997l: 
[ Lebzelter et alj 120001: iLebzelter et al.l [2005: Le bzelter & Wood! 
1201 lh . Our stars show velocities faster than this limit. However, 
an explanation could be that these stars have radial velocities of 
ab out 20 km s" 1 on top of the ir Mira pulsations. 

ILebzelter & Wood! (1201 lh recently investigated the long pe- 
riod variables of NGC 2808. Unfortunately, they had to exclude 
the inner regions of the cluster due to crowding. We identify their 
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long period variable stars in our pointing and found two matches. 
The first one is LW 15 with a period of 332 days while the second 
one, LW 7, does not have a precise position and is therefore an 
uncertain candidate with a period of 21 days. We use the matches 
in our pointing to compare these spectra with our two fast stars. 
The spectra of the fast stars are shown in Figure |5] No conspicu- 
ous features were found in the Mira spectra nor in the spectra of 
the fast stars. The lack of studies of Mira spectra in the optical 
range, especially in the Calcium Triplet region, prevents us from 
making comparisons with other measurements. 

Another argument against the pulsation velocities is the fact 
that optical lines usually form in the upper atmospheres where 
no shock waves can be observed. Also, these optical velocities 
are nearly always di rected inward relative to the stellar center- 
of-mass as found by Woodf ( fl979b . 

It is worth mentioning that we found no asymmetry in any 
of the lines used to measure the radial velocity. Assuming that 
the high velocity would come from an atmospheric effect, such 
as winds or expanding shells, one would expect the lines to be 
either asymmetric like p-Cygni profiles or even split into two 
parts which is often observ ed in Mira variables (eg. lSchorzll992t 
lLamers & Cassinellilll999h . 

We conclude that thermal pulsations are an explanation for 
the high velocities of our two stars seem rather unlikely. The 
range of velocities reached by these pulsations is not high 
enough to explain the peculiar motions of the stars and unlikely 
to be strongly represented in optical lines. 

4. The bound case - Maxwellian distribution 

Having established that the most likely scenario is that the two 
stars do indeed have velocities exceeding more than three times 
the velocity dispersion of the cluster, we investigate possible ac- 
celeration mechanisms of dynamical nature. We first consider 
the case in which we assume that the stars are bound to the clus- 
ter. Under this assumption, the stars could have been accelerated 
either through encounters as giants or in their main-sequence 
stage. For the latter the acceleration can occur through either 
many two-body encounters or a single encounter with a more 
massive object. In this Section we discuss the first scenario and 
the likelihood of observing such an event. 

4.1. Acceleration through uncorrelated two-body encounters 

The first possibility we consider is that the stars acquired their 
high velocities through relaxation, i.e. a sequence of uncor- 
related, distant encounters with other cluster stars over long 
timescales. In such a case they would be bound to the cluster 
and belong to the tail of the three-dimensional Maxwellian ve- 
locity distribution which is established by relaxation. However, 
stars with velocities ~ 3.5<x are quite rare. The chance of detec- 
tion decreases even more if one also considers the probability 
of detecting the high-velocity star in an integral-field unit (IFU) 
since the star has to be very bright to be detected. The final prob- 
ability is composed of the probability of a Maxwellian velocity 
distribution with a certain velocity dispersion producing a high 
velocity star, the number of stars in our pointing and the proba- 
bility of detecting such a star in our pointing. 

In order to make a quantitative statement we calculate these 
probabilities in the foll owing way. We us e the v elocity disper- 
sion profile obtained in Liitzgend orf et alJ ((2012) to estimate <x 
at the position of the two fast stars. From this, we find that the 
stars have velocities of 3.2 and 3.1 times their local velocity dis- 
persion. The probability of finding a star with a total velocity V 



or higher from a three-dimensional Maxwellian velocity distri- 
bution is then given by: 



P(x) = 1 - 



(4) 



with x = V/cr. 

With this information we can calculate how many stars we 
would expect among all the stars in the pointing using the stars 
which have been found in the HST image. In order to obtain the 
detection probability we run Monte Carlo simulations for 10 6 
spaxels. For each simulated spaxel we draw 10 stars from the 
luminosity function. For 9 of them, we assign random total ve- 
locities from a Maxwellian velocity distribution with cr = 13.4 
kms -1 as well as a random spatial orientation. The tenth star 
is always assigned a velocity value of 44 kms -1 . We build the 
spaxel in the following way: for each star we take the template 
spectrum, scale it with the flux of the star, and shift it with 
the amount of its projected radial velocity. Then we combine 
all 10 spectra to a single spaxel spectrum. The velocity of each 
spaxel is measured using the penalized pixel-fittin g (pPXF) pro- 
gr am developed by Cappellari & Emsellem (2004) as described 
in iLtitzgendo rf et al. I (120121) . Figure [6] shows the result of the 
Monte Carlo simulations. Plotted are all the 10 6 measured spaxel 
velocities as a function of the magnitude of the star with V = 44 
km s _1 in the form of a density contour. This shows that the prob- 
ability of detecting a high velocity star strongly decreases with 
the magnitude of the star. 

In order to calculate the detection probability, we set the de- 
tection limit equal to the minimal velocities within the error bars 
(2 kms" 1 ), 42 kms" 1 for FS1 and 35 kms" 1 for FS2. Spaxel 
with velocities below this limit are not considered to be a detec- 
tion anymore. The probability is obtained by deriving the frac- 
tion of stars which are faster than these limits in overlapping 
bins of 100 stars each. The right panel of Figure |6] shows the 
smoothed curve of the probability for FS 1 (magenta) and FS2 
(green). The plot shows that the probability to detect a star with 
a radial velocity of 37 kms -1 in the IFU is higher by almost a 
factor of 5 than a star with 44 kms -1 in the radial component 
and that for stars fainter than my = 17 both velocities are not 
detectable anymore. 

To derive the expected number of stars with this configura- 
tion we multiply the luminosity function of stars in our pointing 
with the probability from the Maxwellian distribution and the 
detection probability from the Monte Carlos simulations and in- 
tegrate over all magnitudes for both stars. As a final result we 
derive: 



AV,fsi = 0.02 stars 
N r ,ps2 - 0.30 stars 



(5) 
(6) 



For further analysis we also derive the pure detection prob- 
ability fi,, which accounts for (1) the geometric probability that 
the velocity component along the line of sight exceeds a certain 
threshold velocity, and (2) the probability that the star is massive 
and bright enough to be detected in the IFU. The probabilities are 
derived in the same way as described above but without multi- 
plying by the probability from the Maxwellian distribution. Thus 
we integrate only the probability function from the Monte Carlo 
simulations multiplied by the number of stars over the entire 
magnitude range. The derived detection probabilities for stars 
like FS1 and FS2 are f bfS1 = 4 x 10~ 4 and / fe , FS2 = 4 X 10~ 3 , 
respectively. 
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Fig. 6. Monte Carlo simulations of recovering the radial velocities of v= 44 km s _1 and v= 37 km s _1 as a function of the magnitude 
of the star. The left panel shows the measured velocities of every spaxel from the 10 6 realizations as a density contour plot. The 
magenta and green lines mark the detection limits of 42 km s _1 and 35 km s _1 for FS 1 and FS2 respectively. A spaxel with a velocity 
lower than 2 km s _1 (which corresponds to the error bars) minus the actual radial velocity is not considered to be a detection anymore. 
The right panel visualizes the probability of detecting a star with a radial velocity of 44 kms _I (magenta) and 37 kms -1 (green) as 
well as their smoothed curves. 



We conclude that neither FS 1 nor FS2 is likely to originate 
from the Maxwellian velocity distribution. The probability for 
a star such as FS2 is higher than that for FS1 by an order of 
magnitude but still too low to explain the observations with this 
mechanism. It is therefore unlikely that FS1 and FS2 acquired 
their high velocities through relaxation and we study other ac- 
celeration scenarios in the following Sections. 



5. Dynamical encounters 

The numbers we derive from the analysis of the Maxwellian dis- 
tribution and observational effects are too small to explain the 
observed stars. Therefore we test for alternative acceleration sce- 
narios due to dynamical three-body encounters of two different 
types of stars (main-sequence stars and giants) as described in 
this Section. 



5.1. Acceleration in the main-sequence stage 

Here we study a scenario in which FS1 and/or FS2 were 
accelerated by a single dynamical encounter at some 
time during the main-sequence stage. To test this pos- 
sibility we perfo rm numerical th re e-body scattering ex- 
perim en ts (e.g. iHills & Fullertonl Il980t iHut & BahcaijH 
1983t ISigurdsson & Phinnevl 119931; iHeggie et al.1 11 996; 
Gvaramadze et al. 2008, 20091; iGvara madze & Gualandris 
20111) involving main-sequence stellar binaries and single 



stars/compact objects of different kinds. The simulations are 
carried out with the sig ma3 package included in the STARLAB 
software environment (McM illan & Hutll 19961) . 

We adopt a m ass of m — 0.8 M H an d a radius of R = 
O.8(M/M ) O 7 R (Hab ets&Heintzelll981h for both stars, and 
we set the characteristic ejection velocity to V e j = 45 km s . We 
consider four types of binary-single star encounters, which are 
detailed in Table [2] (1) encounters with a main-sequence star of 
the same mass, (2) encounters with a 0.8 M white-dwarf, (3) 
encounters with a 1 .4 M Q neutron star and (4) encounters with a 
1OM black hole. We consider a fifth type of encounters, be- 



Table 2. Dynamical encounters between a main-sequence stellar 
binary and a single star. Here MS indicates a main-sequence star, 
WD a white dwarf, NS a neutron star and BH a stellar mass black 
hole. Outcomes of interest for the ejection of one of the MS stars 
include flybys (F), exchanges (E) and ionizations (I). 



Name 


Encounter 


Outcomes 


El 


(MS, MS) + MS 


F, E, I 


E2 


(MS, MS) + WD 


E, I 


E3 


(MS, MS) + NS 


E, I 


E4 


(MS, MS) + BH 


E, I 


E5 


(BH, MS) + MS 


F, E, I 



tween a binary containing a main-sequence star and a stellar 
mass black hole and single stars. Such binaries originate from 
the encounters of stellar binaries with single black holes. 

We first discuss the simulations involving stellar binaries. As 
is typical in three-body scatterings, we use a Monte Carlo ap- 
proach to randomly generate the angles that define the spatial 
orientation of the binary with respect to the incoming star. The 
bina ry eccentricity is randomly drawn from a thermal distribu- 
tion (Heggia [l975l) . having set a maximum value to guarantee 
that the binary components do not come into contact at the first 
pericenter passage. The impact parameter b is randomized ac- 
cording to an equal probability distribution for b 2 with a maxi- 
mum value th at is determined auto matically for each set of ex- 
periments (see lGualandris et al.l2004t for a detailed description). 
The relative velocity between the single star and the binary cen- 
ter of mass is set equal to the cluster velocity dispersion. We vary 
the initial semi-major axis of the binary from a minimum value 
flmin = 0.05 AU set by the radii of the stars to a maximum value 
flmax = 10 AU. 

While stars are treated as point particles in the integrations, 
all inter-particle distances are monitored to identify physical col- 
lisions when the distance between any two stars becomes equal 
to the sum of their radii. Since we are interested in the ejection 
of stars with high velocity, we neglect all encounters that result 
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Fig. 7. Branching ratios versus initial binary separation for the outcomes of encounters between a stellar binary and a main-sequence 
star El (left) or a stellar mass black hole E4 (right). The label "runaway" refers to encounters that produce an escaper with velocity 
larger than 45 km s _1 . 



in a physical collision. While it is a priori possible to eject col- 
lision products during dynamical encounters, the fact that stars 
FS 1 and FS2 lie precisely on the color-magnitude diagram of the 
cluster strongly argues against such a scenario. 

Possible outcomes of three-body encounters include: fiybys 
of the single star past the binary (F), exchanges of the incom- 
ing star into the binary (E), with ejection of one of the binary 
components, ionization (I) of the binary where three single stars 
get ejected, and collisions. For each type of encounter listed in 
Table|2] we consider all outcomes resulting in the ejection of (at 
least) one main-sequence star, and record the ejection velocity. 
Ionizations are rare and, for our set of parameters, only occur in 
the case of encounters with a black hole. In all the other cases, 
the binaries are too hard to be broken apart by the incoming star 
(see Figure |7J. In the case of encounters between equal mass 
stars, the binary is always too hard to be unbound by the incom- 
ing star. 

Figure Q gives the branching ratios, i.e. the fraction of all 
encounters resulting in a particular outcome, as a function of 
the initial binary separation, for model El and E4. The relative 
importance of fiybys, exchanges and mergers is a strong function 
of the binary semi-major axis, with mergers occurring in about 
40% of the cases for the tightest binaries. 

The average velocity of all main-sequence escapers is shown 
in Fig. [8] for all types of encounters. Encounters with a main- 
sequence star produce average recoil velocities significantly 
lower than required for the observed runaways, for all consid- 
ered values of the binary semi-major axis. Encounters with white 
dwarfs and neutron stars produce only slightly larger velocities, 
and reach the required 45kms~ 1 only for the tightest binaries. 
Encounters with a stellar mass black hole, on the other hand, re- 
sult in typical ejection velocities in excess of 45kms~ I for all 
values of a < 0.4 AU. Given that these are average velocities, 
large velocities can occasionally be achieved even for larger sep- 
arations. 

If encounters with stellar-mass black holes lead to ex- 
changes, they are expected to produce a star-black hole binary 
that will then interact with the dominant population of single, 
main-sequence stars. We find that such binaries typically have 
large eccentricities. For this reason, when simulating encounters 
of type E5, we generate binaries with a suprathermal eccentricity 
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Fig. 8. Average ejection velocities of all main-sequence escapers 
for the different types of encounters considered. For model E5, 
we adopt a black hole mass of 30 M Q . The dotted horizontal line 
marks the characteristic velocity of the runaways. 

distribution. The branching ratios for this process are shown in 
Figure [9] Most encounters of this kind result in simple fiybys or 
mergers. The average velocity of escaping main-sequence stars 
is given in Fig.[8]for the case of a 30 M black hole. Only bina- 
ries with separation < 0.1 AU produce average velocities larger 
than 45kms~'. On average, higher black holes masses do not 
result in larger ejection velocities, because a larger fraction of 
encounters results in a collision. 

The rate of ejection of runaway stars like FS1 and FS2 by 
encounters with a compact object can be estimated as 

n = nXV, (7) 

where n is number density of stars in the core, S is the cross- 
section for the process and V is the relative velocity between the 
binary and the black hole. The cross section can be derived from 
the scattering experiments as 

£ = ^Lx/- (8) 
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Fig. 9. Branching ratios versus initial binary separation for the 
outcomes of encounters of type E5, in the case of a black hole of 
30 M . 



represents the maximum impact parameter for the 
process under consideration and f r is the fraction of encoun- 
ters resulting in escapers with velocity larger than 45 km s _1 (see 
Fig.lTTb. The rate for interactions is therefore 



n~ 1.69 xl(T 9 — 



.V 



yr \10 5 pc- 3 /\10AU/ \10kms 



V 



(9) 



Finally, we multiply the event rate % as a function of the 
binary semi-major axis by the distribution function o f the semi- 
major axis taken from Du quennov & Mavorl (119911) . The final 
rates are shown in Figure [10] where we assume a relative veloc- 
ity equal to the dispersion velocity in the core V = 13kms~', 
and a stellar density n = 10 5 pc~ 3 . We obtain the central den- 
sity by extrapolating our density profile (n(0) = 10 6 pc~ 3 ) and 
multiplying this by a binary fraction of 10%. We caution that 
this is a mere working estimate derived from the total estimated 
stellar mass and the structural parameters of the cluster, and % 
scales linearly with the stellar density. The final rate is derived 
by integrating over all a. 

The number of events producing runaway stars with velocity 
> 45 km s -1 is then 
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Fig. 10. Event rates of the different encounters as a function 
of the separation multiplied by the distribution function of the 
semi-major axis. From this the final event rate can be achieved 
by integrating over the values for the semi-major axis. 



In order to estimate an appropriate number of black holes in 
the core of the globular cluster we use results from our own 
N-body simulations which are described in Liitzgendorf et al. 
(2012, in preparation). Our results show that for a model with 
a retention fraction of f ret = 0.3 and N = 128 000 particles 
there are only one or two black holes remaining after 12 Gyr. 
For higher retention fractions such as / let = 1.0 this number in- 
creases to 20 black holes after 12 Gyr. With N - 50 000 stars 
remaining, these values correspond to black-hole fractions of 
fsH = Nbh/N* = 2 x 10~ 5 and 4 x 10~ 4 respectively, equivalent 
to 20 to 500 black holes in a cluster with 1 x 10 6 stars. These 
values are lower limits though, since we are not only interested 
in the number of black holes at the current time, but more so in 
the average number of black holes over the entire lifetime of the 
cluster. Therefore, an average black hole number of Nbh ~ 100 

seems adequate as a first-order estimate. 

Adopting a relaxation time T r i x = 10 8 yr (lHarrislll996l) . we 
find the expected number of runaway stars produced by encoun- 
ters between a main-sequence binary and a stellar mass black 
hole to be: 



Nr^RNnaTtofi, (10) 

where A^bh is the number of BHs to be found in the core, T r \ x is 
the relaxation time of the core (the high velocity of the stars 
should decrease due to two-body relaxation within this time 
scale) and fi, the detection probability as obtained in Section l4~T1 
Estimating the number of black holes in the cluster core is 
not straightforward. There are two factors which play a major 
role. The first one is the black hole retention fraction / let , i.e. 
the fraction of black holes that does not escape the cluster due 
to kick velocities acquired during their formation. The second 
factor is how efficiently equal-mass black holes eject themselves 
from the cluster core by building a subcore through mass seg- 
regation. This pro cess would leave only a few (0-3) black holes 
in the cluster (e .g. Sigurdsson & Hernq uist 1993; Kulkam i et alj 
19931: iMiller & Hamiltor l20"02HO'Learv et al.ll2006l) . However. 
Mackev et alJ (120071) find different results. In their N-body sim- 
ulations, larger numbers of black holes are retained because the 
ejection rate of black holes decreases as the subcluster evolves. 



A^fsi = 0.3 stars (11) 
N,js2 = 3.4 stars (12) 

The fact that the expected number of stars like FS 1 is an or- 
der of magnitude lower than that of FS2 comes from the detec- 
tion probabilities. A star like FS 1 is less likely to be observed due 
to its brightness and the fact that its velocity vector has almost 
no tangential component. The expected number of stars like FS2 
is in good agreement with our observation of the four other fast 
stars in NGC 2808. 

Keeping in mind all the uncertainties that enter these esti- 
mates, we consider both of these numbers consistent with the 
observation of one star like FS1 and one like FS2. 

5.2. The unbound case - Acceleration in the giant phase 

The measured velocities for stars FS1 and FS2 are very close 
to the cluster escape velocity. This makes it hard to establish 
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Table 3. Dynamical encounters between a binary and a single 
star. FS1 and FS2 indicate the giant stars under consideration, 
MS indicates a main-sequence star, and BH a stellar mass black 
hole. Outcomes of interest for the ejection of one of the giant 
stars include fiybys (F), exchanges (E) and ionizations (I). 



Name 


Encounter 


Outcomes 


Gl 


(FS1,FS1) + FS1 


F, E, I 


G2 


(FS2, FS2) + FS2 


F, E, I 


Ml 


(FS1,FS1) + MS 


E, I 


M2 


(FS2, FS2) + MS 


E,I 


M3 


(FS1,MS) + MS 


E,I 


M4 


(FS2, MS) + MS 


E,I 


Bl 


(FS1,FS1) + BH 


E, I 


B2 


(FS2, FS2) + BH 


E, I 



whether the two stars are bound or unbound to the cluster. Here 
we consider the case in which FS 1 and/or FS2 received a recoil 
velocity during a recent dynamical encounter large enough to 
unbind them from the cluster. 

We test the possibility that FS 1 and FS2 were ejected from 
the cluster core by a dynamical encounter by performing numer- 
ical three-body scattering experiments similar to the ones de- 
scribed in Section 15.11 In this case, considering an ejection ve- 
locity of 45 km s _1 and the current projected distance of ~ 0.3 pc 
from the cluster center, we can assume that the encounters that 
accelerated the two stars took place when the stars were already 
in the giant phase. This is an important assumption, and the one 
that distinguishes the encounters described in this Section from 
those in 15.11 Collisions in close encounters affect giant stars 
much more than they affect main-sequence stars due to the larger 
radii, and therefore larger cross-sections, of the giants. 

We consider three main types of binary-single star encoun- 
ters, which are detailed in Table[3] (1) encounters involving only 
giant stars (models Gl, G2), (2) encounters involving giant and 
main-sequence stars (models Ml, M2, M3, M4), and (3) encoun- 
ters between a binary of giant stars and a black hole (BH) (of 
both the stellar and intermediate-mass kind). For each type of 
encounter listed in Table[3] we consider all outcomes resulting 
in the ejection of (at least) one giant star, and record the ejection 
velocity. Since we are interested in the ejection of giant stars 
with large velocities, we neglect all cases of collisions. 

Encounters of star FS1 and FS2 with other giants (case Gl 
and G2), or main-sequence stars (case Ml, M2, M3, M4) result 
in ejection velocities always much lower than 45 km s , for any 
value of the initial binary separation. It thus ruled out that any 
of the two high velocity giants were ejected by encounters with 
stars of the same mass, whether giants or main-sequence stars. 

We also consider the case of encounters with black holes 
(BHs) of different masses. Figure QT| gives the branching ra- 
tios, as a function of the initial binary separation, for model B2. 
The relative importance of mergers and exchanges depends sen- 
sitively on the binary semi-major axis, with mergers occurring 
in more than 50% of the cases for a < 0.2 AU. 

Average ejection velocities as a function of the initial binary 
semi-major axis are shown in Fig. [12] for models B 1 and B2, and 
for different values of the BH mass M BH = 10,20,50, 100 M . 
Velocity distributions for case Bl are shown in Fig. [l3]for BH 
masses of 10,50, 1OOM and the smallest value of the binary 
semi-major axis a — 1 AU. We find that stellar-mass BHs are suf- 
ficient to accelerate star FS2 to velocities of ~ 45 km s for any 
value of the binary separation, and no BH of intermediate mass 
is required in the cluster to explain its velocity. Average ejec- 
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Fig. 11. Branching ratios versus initial binary separation for the 
outcomes of encounters between a stellar binary and a black hole 
of 10 M (model B2). The label "runaway" refers to encounters 
that produce an escaper with velocity larger than 45 km s _1 . 
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Fig. 13. Distribution of ejection velocities in encounters between 
a binary of giant stars and a black hole (case Bl). The three 
panels refer to BHs of different mass: 10,50, 100 M , from left 
to right. The binary separation is set to 1 AU. 



tion velocities for star FS 1 are somewhat lower, due to the much 
larger size and therefore to the much larger number of close en- 
counters that end with a collision. However, if we look at the dis- 
tributions of ejection velocities for model Bl, which are shown 
in Fig.[l3]for the case a — 1 AU, we find that a 10 M black hole 
is sufficient to produce ejection velocities of 40 - 45 km s _1 . 

We compute the ejection rates for this ejection mechanism 
in the same way as in Section IBTTI To obtain the final number of 
ejected stars, however, the relaxation time in Equation [8] needs 
to be replaced with the crossing time of the core. With a typical 
velocity of 13kms _1 , the crossing time for a core diameter of 
1 .4 pc is T CI ~ 10 5 yr. This results in a final number of events as 
listed in Table |4] It is clear that the higher the black hole mass, 
the more stars with velocities larger than 45 km s _1 are produced. 
However, even a black hole of 10 3 M Q would not be enough to 
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Fig. 12. Average ejection velocity in runs Bl (left) and B2 (right) as a function of the initial binary semi-major axis. Different lines 
refer to different values of the BH mass, Mbh = 10,20,50, 100 M , respectively, from bottom to top. Dotted lines indicate the 
required ejection velocity of 45 km s _1 . 



Table 4. Expected numbers of stars if accelerated by an en- 
counter of a giant and a massive black hole. 







NUMBER OF EXPECTED STARS 




10 M 


20 M 50 M Q 100 M e 


FS1 


3.17 x 10 


6 1.46 x 10~ 5 3.38 x 10~ 5 3.77 X 10~ 4 


FS2 


1.70 x 10 


" 5 2.21 x 10~ 5 3.04 x 10~ 5 6.37 x 10~ 5 



produce an acceptable number of stars to explain the observa- 
tions. 

We conclude that a dynamical encounter between a binary 
of giant stars and a stellar-mass black hole of > 1OM can 
produce giant stars with velocities ~ 45kms~'. However, the 
number of expected fast stars produced by this mechanism is one 
to three orders of magnitudes lower than our observations show. 
The reasons for this discrepancy are, on the one hand, the lower 
ejection rates due to the shorter timescales in which the stars 
will disappear from the field of view and, on the other hand, the 
fact that giant stars are less likely to achieve large velocities than 
main-sequence stars. For these reasons, this scenario is unlikely 
to have accelerated the observed runaways. 

6. Comparison with other globular clusters 

In Section[T]we me ntioned the cases of M3 and 47 Tuc and their 
high-velocity stars. Meylan et al. (1991) argued that the two fast 
stars in 47 Tuc, which move with 3.6cr c and 4.8<x c , must have 
been recently accelerated by encounters with a giant star, sim- 
ilar to our unbound case. The reason for this assumption is the 
fact that the velocities of these stars (Vi = -36.7 kms -1 , V2 = 
32.4 km s" 1 ) are, wi thin the errors, equal to the escape velocity 
V esc = 35.38 km s _1 dPeterson & Kingf l975). This is a very sim- 
ilar situation to our case and therefore we would be more cau- 
ti ous in claiming the st ars are bound or unbound. The argument 
of Meyla n et al. d!991l) that these stars could only recently have 
been accelerated is questionable. Furthermore, even if the veloc- 
ities of the stars exceed the escape velocity by a few km s , the 
probability that the stars are able to actually escape the cluster 
immediately remains low, since they are only able to escape by 



passing through small apertures near the Lagrange points LI and 
L2. This can result in a significant delay of the escape of stars 
dFukushige & H eggie 2000). For this reason it is a valid hypoth- 
esis that these stars also got accelerated in their main-sequence 
stage. Because the central density[] k>g(po/ pc~ 3 ) = 4.88 and core 
relaxation time log(f e / yr) = 7.85 (lHarrislll996l) of 47 Tuc are 
very similar to the values of NGC 2808, the expected numbers 
of fast stars will be close to our derived numbers. 

In the case of M3 the central velocity dispersion is lower than 
the one of NGC 2808 by 50 % (cr c = 4.9 kms" 1 ). However, it 
is important to note that the absolute velocity dispersion is not 
the key factor, but r ather its value relative t o the velocities of the 
high-velocity stars. iGunn & Griffinl (fT979h measured velocities 
of 17.0 kms -1 and -22.9 kms -1 for two stars which are most 
likely cluster members. Again, these velocities lie in the same 
range as the escape velocity (V esc = 20.81 kms -1 ) but, as in 
the two previous cases, do not exceed this limit by much. The 
only difference with NGC 2808 and 47 Tuc is the central den- 
sity of log(po) = 3.57. This would lower othe expected num- 
ber of high-velocity stars by an order of magnitude. Since this 
would most strongly affect the prediction for the unbound case, 
the most likely scenario for M3 is still that the stars are bound 
and were accelerated by a dynamical encounter with a massive 
object while they were in the main sequence phase. 

The globular cluster NGC 6752 shows two millisecond pul- 
sars outside the core: a recycled binary pulsar at more than three 
half-mass radii and a single pulsar at 1.4 half-mass radii. Such 
systems sugge st the presence of a mas sive object in the center 
of the cluster. IColpi et alj (120021 120031) show that the off-center 
location of the pulsars can be explained with scattering from 
a binary of stellar mass black holes or a single intermediate- 
mass black hole. This supports the idea that high-velocity stars 
in globular clusters may be produced by single or binary black 
holes. 



1 With the caveat that log(po) is measured with a different technique 
than our derived value and therefore lower than our actual central den- 
sity for NGC 2808. 
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7. Summary and conclusions 

We detect five high-velocity stars in the globular cluster NGC 
2808 using integral field unit spectroscopy as well as Fabry- 
Perot data for individual stars. All stars are ~ 0.3 - 0.5 pc 
away from the cluster center and their velocities correspond to 
3.6 - 4.1cr c , where cr c is the velocity dispersion in the core. All 
stars are very likely to be cluster members due to their positions 
in the CMD. Furthermore, we find that the proper motions of the 
stars are in good agreement with the proper-motion distribution. 

For the case of the two fast stars found in the integral-field 
unit, we discuss various scenarios which could explain their pe- 
culiar motions. By performing a binary synthesis and looking at 
HST far UV and near UV images, we exclude the possibility of 
these stars being close binaries. Furthermore, we discuss the pos- 
sibility of atmospherically acitve stars such as long period vari- 
ables and conclude that this scenario can also be excluded due to 
a lack of spectroscopic evidence and general inconsistency with 
the observed velocity range. 

The measured velocities of the stars are close to the esti- 
mated escape velocity from NGC 2808. For this reason, we con- 
sider two cases of acceleration: The bound case, in which the 
stars were accelerated while in the main-sequence stage; and 
the unbound case, in which the acceleration must have taken 
place recently. We perform numerical three-body scattering ex- 
periments for both cases. In the bound case we also test the pos- 
sibility that these stars belong to the tail of the Maxwellian ve- 
locity distribution by means of Monte Carlo simulations. In the 
case of three-body encounters of main-sequence stars and com- 
pact objects (the bound case) we consider four types of encoun- 
ters: (1) encounters with a main-sequence star of the same mass, 
(2) encounters with a 0.8 M white-dwarf, (3) encounters with 
a 1.4M neutron star and (4) encounters with a ~ 10M G black 
hole. To obtain the expected numbers of ejected stars, we cal- 
culate the event rates of each encounter producing a star with a 
velocity larger than 45 kms~' as a function of the binary sepa- 
ration. Combining the r ates with the log-normal orb ital period 
distribution obtained by Duquennov & Mavorl (ll99ll) . we obtain 
the total event rate by integration. 

Similarly, we test for the possibility of encounters with gi- 
ant stars (the unbound case). As an input, we take the estimated 
properties of the two stars (radius and mass) and consider three 
main types of binary-single star encounters: encounters between 
giants, encounters of giant and main-sequence stars, and encoun- 
ters including a black hole. We find that all encounters including 
only stars result in much lower escape velocities than the ob- 
served ones. Encounters between a binary of giants and a black 
hole, however, result in velocity distributions which peak at the 
observed velocities. Especially a black hole of Mbh ~ 10 M 
and a binary separation of a = 1 AU predicts velocities around 
40 kms -1 for star FS1. For star FS2, any separation produces 
large enough ejection velocities. 

The scattering experiments show that no intermediate-mass 
black hole is needed in order to accelerate the stars to their veloc- 
ities. Furthermore, only the encounters of main-sequence stars 
and compact objects have rates high enough to explain the ob- 
served stars. This leads us to the conclusion that the stars are 
bound and were accelerated by an encounter with a ~ 1OM 
stellar black hole. This cluster must thus have a high concentra- 
tion of stellar-mass black holes in its center in order to make 
such an encounter likely enough to be observed. This is in 
agreement with the latest results on the central kinematics of 
NGC 2808 where no int ermediate-mass black hole is detected 
Liitzgend orf et alJ (|2012). Such a retention of stellar remnants in 



evolving globular clusters is also in line with theoretical expec- 
tations, which allow up t o 30% - 50% of the cluster mass to be 
const ituted by remnants (Baumgardt & Makino 2003; Kruiissen 
I2009h . 

We compare our results with the two other known cases in 
which high velocity stars were detected in globular clusters, 
M3 and 47 Tuc. These cases are very similar to NGC 2808 
and might therefore be explained with the same arguments. The 
high-velocity stars in these clusters are most probably products 
of encounters of main-sequence binaries with a massive black 
hole. For further analysis, N-body simulations could help veri- 
fying these results by running specific simulations dedicated to 
the statistics of fast stars and their origin. When combined with 
an analysis as carried out in this paper, this provides a general 
prediction for the origin of high-velocity stars in globular clus- 
ters. 
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